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Abstract 

The physics of a two-level system coupled to Ohmic bath is studied by means of the 
perturbation approach based on a unitary transformation. Our main results are: The 
coherence- incoherence transition is at a c = |[1 + A r /u c ]; for a < a c the dynamical 
quantity P(t) = cos(uj t) exp(— jt); the sussceptibility x"(u;)/u; is of a double peak 
structure for a < a c and the Shiba's relation is exactly satisfied; at the transition 
point a = a c the real time correlation function C(t) w —l/^t 2 in the long time limit. 
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Studies of the two-level system coupled to Ohmic bath (spin-boson model, SBM) 
have a long history]]], 0. Although there were many studies using various kinds of 
methods, in recent years the physics of SBM has still attracted considerable attention 
because it provides a universal model for many physical systems]]], The Hamiltonian 
of SBM reads 

H = --Aa x + ^u k b\b k + -^g k {b\ + b k )a z . (1) 

~ k k 

The notations are the same as usual]]], 0. In this work we consider the zero bias case 
with temperature T = 0. The Ohmic bath is characterized by its spectral density: 
J^kQk^i^ ~ u k) = 2au9(u c — u>), where a is the dimensionless coupling constant and 
6{x) is the usual step function. 

The Hamiltonian (1) seems to be quite simple. However, it cannot be solved ex- 
actly and various approximate analytical and numerical methods have been used[l-13]. 
Recent studies of SBM focused on its dynamical properties[6-12]. Among these stud- 
ies, the variational approach proposed by Silbey and Harris |H| is of interest. They 
proposed to make a unitary transformation: H' = exp(S)H exp(— S), with 

s = Y.frUA-h)cr z . (2) 

k k 

Then, the variational ground state energy is: E g = (si\({0 k }\H'\si)\{0 k }) , where \si) = 
^ | ^ and |s 2 ) = ( ^ ^ are eigenstates of a x and |{0fc}) is the vacuum state 
of bosons in which n k = for every k. By minimizing E g £ k is determined as 

6; = uj k /(u k + 7) A). (3) 



The authors calculated the renormalized tunnelling A r ]T3" 



A r ( A\ — 



V=^=[e-J exp 



a ( r]A u c + r]A" 

+ In ■ 



(4) 



1 — a \ uj c + rj A uo c 

Although the variational method can predict correctly the delocalized- localized tran- 
sition point ati = 1 at the scaling limit A/u c 1, it is not suitable for calculating 
the dynamical properties of SBM. In this work we present a new analytical approach, 
based on the unitary transformation ((2) and (3)) and the perturbation theory [Tjj], for 
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calculating the dynamical preperties of SBM which works well for the coupling con- 
stant < a < 1 and the bare tunnelling < A < u c . The approach is quite simple and 
physically clear, and may be easily extended to more complicated coupling systems. 
Throughout this work we set h — 1 and ks — 1- 

The unitary transformation (2) can be done to the end and the result is H' = 
Hq + H[ + H 2 , where 

H' = -± V Aa x + J2io k bth - E £^( 2 - (5) 

k k k 

^=«p[-E^a. ( 6 ) 
k iuj k 

H[ = \Y, 9k(l - tk)(b\ + h)a z - ^rjAiay E ~Ub\ - h), (7) 
H' 2 = - l -Aa x fcosh{E ~UA ~ h)} ~ v) 
-Uia y fsinh{E 9 -Ub\ - b k )} - vJ2 -UA ~ bk)) ■ 



(8) 



Obviously, H' Q can be solved exactly because in which the spin and bosons are decou- 
pled. The eigenstate of H' is a direct product, |s)|{nfc}), where \s) is \si) or |s 2 ) and 
|{rifc}) means that there are n k phonons for mode k. The ground state of H' is 

lflb> = |si>|{0*}>. (9) 

H[ and H 2 are treated as perturbation and they should be as small as possible. Eq.(3), 
ik = u k /(u k + r]A), leads to 



H i = 7^ A E ~ & b\{a z - io y ) + b k (a z + ia y ) 
* k uj k l 



(10) 



and H[\g ) = 0. Note that the form of Eq.(3), in this work is determined by 
H[\g ) = 0, instead of minimizing the ground state energy E g . Thus, we get rid of the 
variational condition. 

The lowest excited states are | s 2 ) | {0 fc }) and where \l k ) is the number state 

with n k = 1 but n k > = for all k' ^ k. It's easily to check that (go\H 2 \go) = 
(because of the form of rj in Eq.(6)), ({0^)1 (s 2 \H 2 \g ) = 0, (lkKs^H^lgo) = 0, and 
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{{Qk}\(s2\H' 2 \si)\l k ) = 0- Moreover, since H[\g ) = 0, we have ({0 k }\{s 2 \H[ \g ) = 
and (si\H[\g ) = 0. Thus, we can diagonalize the lowest excited states of H' as 

H' = —-i]A\g ) (g \ + ^2 E\E) (E\ + terms with higher excited states. (11) 

2 E 

The diagonalization is through the following transformation J3J: 

h>|{o fc }> = £*(£)|£>, (12) 

E 

\s l )\l k ) = J2yk(E)\E), (13) 

E 

\E) = J2x(E)\s 2 )\{0 k }) +J2y k (E)\ Sl )\l k }, (14) 



where 



_ , , ~ 1/2 

' V(t) ■ ' ^ (E + v A/2 - u k )\ 



V 2 



(15) 



**W = TP -L /T/9 ( 16 ^ 

E + 7] A/2 - u k 

with Vfc = rjAg k ^ k /uj k . E's are the diagonalized excitation energy and they are solutions 
of the equation 

E---Y,T-^ = 0- (17) 

2 ^ E + V A/2-co k 1 ; 

The dynamical quantity P(t) = (b, +1| (+l\e lHt a z e~ lHt \ + 1)|6, +1) is defined in 
Ref. 0, where | + 1) is the eigenstate of a z = +1 and |6, +1) is the state of bosons 
adusted to the state of a z = +1. Because of the unitary transformation (e s a z e~ s = o z ) 

P(t) = ({0 fe }|(+l| e ^'V z e-^' f | + l)|{0 fe }), (18) 

since e 5 | + 1)|6,+1) = | + l)|{0 fc }). Using Eqs.(ll)-(17) the result is 

P(t) = ^5> 2 (£)exp[-i(£ + 7?A/2)£] + ^ J> 2 (£) exp[i(£ + r?A/2)t] 

E 2 E 

1 r ( t/2 \ _1 



dE'e-^ [E'-r,A-J2 — 



_k 

Am Jc V" ' _ ^E' + i0 + -uj k ) 

r2 \ "I 



+ -L ( f d E>e iE ' t [E'-r ] A-y] \ ] , (19) 

AmJc V k E -iO + -uj k J ' v ' 



where a change of the variable E' = E + rj A/2 is made. The real and imaginary parts 
of Efc Vk/( E> ± - wjfe) are denoted as J?(£") and T7(£')> 

|w|(w c + r/A) 

Ml 

lu + r/A I uo c + r/A + r/A 



(20) 



r]A(u; c - 

7(V) = 2a7rcj(r7A) 2 /(cj + r/A) 2 . (21) 

The integral in (19) can proceed by calculating the residue of integrand and the result 
is P{t) = cos(u;oi) exp(— jt), where Uq is the solution of equation 

uj - r]A - R(u) = (22) 

and 7 = 7(77 A) = a r nr]A/2 (the second order approximation). This P(t) is of the 
form of damped oscillation and one can check that the solution ujq is real when 1 > 
2au c / (uj c + 77 A). When 1 < 2au c / (uj c + 77 A), the solution u is imaginary and we have 
an incoherent P(t). a c = |[1 + r]A/u> c ] determines the critical point where there is a 
coherent-incoherent transition. Note that when A/u c <C 1, we have a c = 1/2, ujq = 
and P(t) = exp(— j c t) (7 C = 7reA 2 /4u; c since 77 = eA/uo c from Eq.(4)), which is the 
same as was predicted by previous authors (7 C = 7rA 2 /2uj c in Ref.|| £|, §]). Fig.l shows 
the calculated uq/ A r as functions of a (a < a c ) for A = 0.01, 0.1, and 0.5. A r = r/A 
is the renormalized tunnelling. The dotted line is 7/A r = air/2. 
Since e s a z e~ s = a z , the retarded Green's function is 

G(t) = -i9(t) ([exp(iH't)a z exp(-tH't),a z }}' , (23) 

where (...)' means the average with thermodynamic propability exp(—/3H'). The 
Fourier transformation of G(t) is denoted as G(u), which satisfies an infinite chain 
of equation of motion|]15|. We have made the cutoff approximation for the equation 
chain at the second order of g k and the solution at T = is 

1 1 
G ^ ~ uJ-r,A-Y. k V£/{uj - uj k ) ~ u + r,A-Y,kVi/{uJ + u k y (24) 

The susceptibility x( w ) = ~~ G(u), and its imaginary part is 

V i ( A = j{u)8{u) <y(-u)0(-u) 

Xy> [u - r/A - R(u)] 2 + 7 2 H [cu + 77A + J R(-w)] 2 + 7 2 (-^)' y ' 



The uj — > limit of S(u) = x' '(w) / 'uj is 



uj (r?A) 2 {l - 2ct[l -r/A/(cj c + r/A)]} 2 ' ^ 



Besides, the real part of the susceptibility is 



/ 2 
X ^ = 0) = V A{1 - 2a[l - 7] A /{u c + r] A)]}' (2?) 



Thus, the Shiba's relation 0, ft g |2j 



lim^l = W(c. = 0) 2 (28) 
w->o 2 



is exactly satisfied. 5 (a;) has a double peak structure for a < a c . For a > a c there is 
only one peak at u = 0. Fig. 2 shows the 5(o;) versus u; relations for fixed a = 0.3 and 
A = 0.01, 0.1, and 0.5. 

The symmetrized correlation function 

C(t) = ^Tr{exp(-0H)[a z (t)a z + a z a z (t)]}m^M-W} 



1 f°° (3uo 

— I do;coth(— )lmG(u) exp(— tut) 
2n J-oo 2 



1 r 00 

= ~ / duJ l A J M2i 2 COS M)' ( 29 ) 

n Jo [uj — r]A — it(a;o)J + 7 

where ujq is the solution of Eq.(22) and 7(0;) in the denominator is approximated by 
the second order approximation 7. At the scaling limit A/uj c <C 1 and the coherence- 
incoherence transition point a = 1/2, 

c{ t )= r ^T^, ( " A L - (3°) 

w 7o cu 2 + 7 C 2 (a; + ??A) 2 v ; 

C(t) decays algebraically in the long-time limit: C(t) ~ — l/7 2 t 2 , which is the same as 
what was predicted by previous authors. 

In summary: The physics of SBM is studied by means of the perturbation approach 
based on a unitary transformation. Analytical results of the dynamical quantity P(t), 
the sussceptibility x"( u ) an d the real time correlation function C(t) are obtained for 
both the scaling limit A r /uj c <C 1 and the general finite A r /u c case. Our approach is 
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quite simple, but it can reproduce nearly all results which agree with those of previ- 
ous authors using various complicated methods. Besides, our approach can be easily 
extended to other more complicated coupling systems. 
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Figure Captions 



Fig.l w /A r versus a relations for A/u c = 0.01 (solid line, a c = 0.50014), 0.1 
(dashed, a c = 0.51212), and 0.5 (dashed-dotted, a c = 0.66244). The dotted line is 
7/A r = avr/2. 

Fig. 2 S(lu) as functions of uj for fixed a = 0.3 and A = 0.01 (solid line), 0.1 
(dashed), and 0.5 (dashed-dotted). 
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Fig.2 



